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Mass-gap for black hole formation in higher derivative and ghost free gravity 
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We study a spherical gravitational collapse of a small mass in higher derivative and ghost free 
theories of gravity. By boosting a solution of linearized equations for a static point mass in such 
theories we obtain in the Penrose limit the gravitational field of an ultra-relativistic particle. Taking 
a superposition of such solutions we construct a metric of a collapsing null shell in the linearized 
higher derivative and ghost free gravity. The latter allows one to find the gravitational field of a 
thick null shell. By analysing these solutions we demonstrate that in a wide class of the higher 
dimensional theories of gravity as well as for the ghost free gravity there exists a mass gap for the 
mini black hole production. We also found conditions when the curvature invariants remain finite 
at r = 0 for the collapse of the thick null shell. 
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It is generally believed that the theory of general rela¬ 
tivity (GR) should be modified to improve its ultraviolet 
(UV) behavior and remove singularities. One of the op¬ 
tions is to allow terms in the gravitational action that 
contain more than two derivations. The UV properties 
of the higher derivative theory of gravity are usually bet¬ 
ter than in GR. In particular, the forth order gravity 
can be made renormalizable |l[. At the same time, the 
gravitational potential of a point mass in the Newtonian 
limit of such theories is usually finite (see e.g. as and 
references therein). However the higher derivative grav¬ 
ity possesses new unphysical degrees of freedom (ghosts) 
1), ;2| . The problem of ghosts can be solved if one allows 
an infinite number of derivatives in the gravity action, 
that makes it non-local. Ghost-free theories of gravity 
are discussed in M- Their application to the problem 
of singularities in cosmology and black holes can be found 
in S 

In this paper we study gravitational collapse of a small 
mass in higher derivative (HD) and host-free (GF) the¬ 
ories of gravity. We obtain solutions of the linearized 
equations for such theories for a spherical collapse of null 
fluid. We demonstrate that if a static gravitational field 
of a point mass in the HD and GF gravity is regular at 
r = 0 [Ij, then the metric for the collapsing object has 
the same property. This means, that the perturbation of 
the metric, which is proportional to the collapsing mass 
M , is smooth and uniformly bounded, so that the higher 
in M corrections can be neglected in the leading order. 
This implies that for the collapse of a small mass an ap¬ 
parent horizon is not formed. In other words, for this 
wide class of HD and GF theories of gravity there exists 
a mass gap for mini black hole production. This property 
is a consequence of the existence of the UV length scale, 
where such theories become different from GR. For the 
Weyl modified gravity this was shown long time ago in 

i- 

We study the linearized gravity equations on the flat 
Minkowski background tj^ and write the metric in the 
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form = rj^ + h The most general action for the 
higher derivative theory of gravity which contains not 
higher than second power of h is jH, @ 


S 



\ h„ v a(D)Dh^ + h;b{U)d a d v h^ 


+ hdp)d ll d v h lu ' + ^hd(D)Dh 

+ h X<rfM da d X d^d v h^ . 


(1) 


Five, in general non-linear functions of the box-operator 
obey the following three relations 


a + b = 0 , c + d = 0 , b + c + f = 0 . (2) 


Thus the action S contains in fact only two independent 
arbitrary functions of the box operator. In order to re¬ 
cover GR in the infra-red domain these functions must 
satisfy the following conditions a(0) = c(0) = —6(0) = 
—d( 0) = 1. Let us list some special interesting examples 
0 : 

1. General relativity (GR), C = R: a = c = 1; 

2. Gauss-Bonnet (GB) gravity, C = R+a(\3)Q , where 
g = R 2 - AR^R^ + R^ a pR^ afi is the Gauss- 
Bonnet invariant: a = c = 1; 

3. L(R) gravity, C(R) = £(0)+£'(0)R+l/2£"(0)R 2 + 
...: a = 1, c = 1 — £"(□); 

4. Weyl gravity, C = R- a = 1- 

g~ 2 D, c = 1 — 

5. Higher derivative (HD) gravity: a = J([" =1 (l — 
u~ 2 0), c = TIC =1 (! ~ v k* u )- For simplicity, in 
what follows we assume that masses fit are differ¬ 
ent; 

6. Ghost free gravity: a = c = exp(— D/g 2 ). 

It is evident that the linearized GB gravity has the same 
properties as GR. In the same linearized approximation 
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the Weyl and L(R ) theories of gravity are nothing but 
special cases of the general HD gravity. 

Let us consider first static solutions of the linearized 
gravity equations. In the Newtonian limit the stress- 
energy tensor is r = p(f)( 5°<5°, and the metric is of the 
form 


ds 2 = — (1 + 2ip)dt 2 -f- (1 — 2-0 -f- 2<p)d£ 2 . (3) 

The functions tp and ip obey the equations U] 

a(A) A ip = 8 ttGp , (4) 

(a(A) — 3c(A))(Ay> - 2 A i/j) = 8irGp. (5) 

Here A is a usual flat Laplace operator in a flat 3D space 
with metric di 2 , and G is the gravitational coupling con¬ 
stant. After solving (© and finding the potential ip, one 
can find the second potential ip by solving ©• 

For a point mass p = md(r) the solution ([3]) is spher¬ 
ically symmetric. We call it finite if p(r) and ip{r) near 
r = 0 have the form 

ip(r) ~ip 0 + ipir + ^ip 2 r 2 +0(r 3 ), (6) 

<p{r) ~ Po + Vir + + 0(r 3 ). (7) 

A finite solution is not necessary regular one. Really, the 
Kretschmann tensor invariant 1Z 2 = R a p^gR 01 ^ 6 for the 
metric ©, © and © is the form 

n 2 = t! + V + O(i), (8) 

A 2 = 8(4 ipl - bipiipi + 3<fl), 

A\ = 16[V>i (5ip 2 - 4^ 2 ) - 4(pi (ip2 ~ <£ 2 )] • 

The quantity A 2 is a positive definite quadratic form of 
variables ipi and p\, and it vanishes only when ip\ = 
p\ = 0. In such a case the quantity A\ vanishes as well, 
so that 1ZP 2 is finite at r = 0. We call such a solution 
regular. We also call a solution 0-regular, if ipi = 0. For 
a special class of theories, where a = c, one has ip = 2<p 
and a solution which is 0-regular is at the same time a 
regular one. 

We denote O = a(A)A, 

Q(0 = 0 -1 (A = —£) = —[£ci(—£)] _1 • (9) 

and assume that Q(£) can be written as the Laplace 
transform of some function /(s) 

/>00 

Q(0 = / dsf(s)e " sf , (10) 

JO 

1 ncx.-\-ioo 

/(^) = W— / ^Q(£)e s «. (11) 

‘-‘Al J QL — ioQ 


1 Notice that the first of the equations (17) of this paper contains a 
misprint and it should be written as 2(a—3c)[V 2| t > —2V 2 '!'] = rep. 


The second relation is nothing but the inverse Laplace 
transform. A parameter a must be chosen so that the 
integration path in m lies in the domain of the analyt- 
icity of Q(£). 

A formal solution of the operator equation OG = — I 
can be written by using the Laplace transform GB¬ 
it contains the exponent exp(sA), which in the x- 
representation is nothing but the heat kernel 

„ — \x—x'\ 2 / (4s) 

{x'\e sA \x) = K{\x-x'\\s) = ( 47 rg) 3/2 ■ ( i2 ) 


Thus the potential 0(r) for a point mass is 

POO 

ip{r)=8T:Gm / dsf{s)K{r\s) 
Jo 

Gm r a+io ° 

= — / d£Q(£) 

'K'l'f' Ja—ioo 




(13) 

(14) 


We consider at first a case of HD gravity. We assume 
that the function Q(£) has simple poles and write it in 
the form 


Q(£) = -[en( 1+ ^)] _1 - ( 15 ) 


This covers the above listed gravitational theories 1-5 
except for some degenerate cases. 

The Heaviside expansion theorem 0 gives the follow¬ 
ing expression for /(s) 


/oo = -( i-E/r 1 ^ 4 ), (i6) 

i =1 

where Pi = TI™=i^,(1 - Taking the integral 

Cl one obtains 


ip(r) =—2Gmr X (1 — ^ 1 e /iir ), (17) 

i =1 

For GR f(s) = 1 and one has 

1 p{r) = 2 ip(r) = —2Gm/r . (18) 

For a theory with higher derivatives, where n > 1, the 
potential ip(r) near r = 0 has a form © with 

n 

ipo = —2GmSi , iP 1 =GmS 2 , Sfc = ^Ab^” 1 • (19) 

i=l 

We used here a relation So = 1. Thus a theory with 
higher derivatives is '0-regular, when the condition S 2 = 0 
is satisfied. 

For the GF gravity/(s) = — d(s — p~ 2 ) and one repro¬ 
duces the result of 

0(r) = 2ip(r) = —2 Gm er£(pr/2)/r. 

This solution is regular at r = 0. 


( 20 ) 
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We demonstrate now, how using a solution of ((4]) for a 
static point mass one can obtain a solution for an ultra- 
relativistic particle. Let us write the flat space metric in 
the form dt 2 = dy 2 + d( 2 . and suppose that the source, 
generating the gravitational field ([3]), moves along the y- 
axis with a constant velocity /?. To find the gravitational 
field of the moving source we make the following boost 
transformation 

t = \-V + \+u, y = X-V — X+u. (21) 

Here X± = (1 ± ( 8 ) 7/2 and 7 = (1 — /3 2 ) -1 / 2 . In the limit 
7 —> 00 one gets y ~ —7 u, t ~ 7 u, £ 2 ~ 7 2 u 2 + Ci> and 


ds 2 = —dudv + d( 2 + dh 2 , 

dh 2 = &du 2 , $ = —2 lim (y 2 ^). ( 22 ) 

7—^00 

We assume that the energy of the particle, M = 7 m 
remains constant in this (Penrose) limit. We use also the 
following relation 

lim 7 exp(—7 2 m 2 /( 4 s )) = y/4ns8(u ). (23) 

7—^OO 

These relations and m give 

$ = -4 GMF{( 2 ± )6(u) , (24) 

F(z) = _/( s )e- z /( 4s ). (25) 

Jo s 


For GR (as well as for GB and L(R) gravity) one has 
f(s) = 1 and F(z) = ln(z/r] 2 ), where 77 is an infrared cut¬ 
off parameter. The relations (l22ll and (l24l) correctly re¬ 
produce the well known Aichelburg-Sexl solution for the 
gravitational field of an ultra-relativistic particle (’’pho¬ 
ton”) in GR. 

Using expression (fTbl) for f(s) for the HD gravity and 
taking integral in (l25l) one finds 

n 

F(z)=\n(z/r] 2 ) + 2'Y^p- 1 K 0 (ii i y/z). (26) 

i=1 

In the limit /ij —> 00 the second term in the right-hand 
side vanishes and one obtains the correct expression for 
GR. In the presence of the higher derivatives the leading 
term of the function F(z) at small z is 

F(z) ~ C - ^S 2 z(\nz - 2c) - J Sz + 0(z 2 ), (27) 

where c = 1 + In 2 — 7 and S = V) ",, Hi ln^DP^ -1 . For 
the ghost-free gravity one has [ill ] 

F(z) = In 2 + 7 + Ei(l, z) ~ ^ - jz 2 + 0(z 3 ). (28) 

The obtained metric (l22l) . (l25l) can be used to find a 
solution for the linearized HD and GF gravity equations 
for a collapsing spherical thin null shell. For this pur¬ 
pose one considers a set of ’’photons”, passing through a 


fixed point P of the Minkowski spacetime. In the contin¬ 
uous limit this set fills the surface of the null cone, with 
the vertex at P. We additionally assume that the den¬ 
sity of this spherical distribution of the ” photons” is uni¬ 
form and the corresponding mass per a unit solid angle is 
M/ 47 t. Since we are working in the linear approximation 
the resulting gravitational field for such a distribution is 
ds 2 = dsy + (dh 2 ), where (dh 2 ) is obtained by averaging 
of a single ” photon” metric over their spherical distribu¬ 
tion. The calculations give jTl] 

ds 2 = - dt 2 + dr 2 + r 2 duj 2 + (dh 2 ), z = r 2 - t 2 , 
(dh 2 ) = —2GMr~ 1 F(z)[(dt — tdr/r ) 2 + zduj 2 /2] .(29) 

Let us denote 

9 = (Vp ) 2 , p 2 = g es = r 2 - ^-zF(z ), (30) 

then the equation <? = 0 determines a position of the 
apparent horizon, if the latter exists. In the linear in M 
approximation this function is 

g = l-2GMr~ 1 q(z), q(z)=zF\z ), (31) 

where (...)' = d(.. .)/dz. For GR (as well as for GB and 
L(R) gravity) q(z) = 1. 

Using (l26l) one finds that for the HD gravity 

n 

q(z) = 1 - y/zy^JljPj ' 1 K 1 (// jy/z) . (32) 

i =1 

For small z one has 

q(z) = — -^S 2 z(\nz — 2c + 1) — —Sz + 0(z 2 ). (33) 

Let us demonstrate now that the function g is positive 
for small enough M, and, hence, the apparent horizon 
does not exist. Let us notice t hat out side the null shell 
\t\/r < 1. We denote t = ±y/l — /3 2 , 0 < /3 < 1, then 
one has (y t = (dp. t r) 

n 

q(z)/r = p'Y^H i P i 1 Z(y i ) , Z(y) = - - K^y ). (34) 

i -1 y 

The function Z(y) is positive and takes maximal value 
0.399 at y = 1.114. Thus 

n 

W(z)\/r < 0.4^/iilPil" 1 . (35) 

i=1 

This implies that for small enough value of the mass M 
the invariant g is positive everywhere outside the shell. 
In other words, for such mass M the collapse of the null 
shell does not produce a mini black hole. This means 
that for the class of the higher derivative theory of gravity 
(1151) with n > 1 there is a mass gap for the mini-black 
hole production. The value of this map is determined 
by the characteristic length scale p~ 1 of the theory. The 
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apparent horizon does not exist if Mfi < 1. The same 
conclusion is valid for the GF theory of gravity Ha¬ 
lt is possible to calculate the curvature invariants for 
the metric (1291) . In particular, the Kretschmann curva¬ 
ture invariant TZ 1 2 3 in the lowest order in M is 

TZ 2 = 48< ^6 M V , T = 2 z 2 q' 2 - 2 zqq' + q 2 . (36) 

Using expression (l33l) for small 2 one finds 

T ~ ^ 2 [{w 2 +4 w + 5 )S 2 + 2 (w + 2)SS 2 + S 2 }, (37) 

where w = In z — 2c. This means that Kretschmann cur¬ 
vature vanishes on the null shells. However, in a general 
case it is divergent at r = 0. 

We demonstrate now that this divergence is a result 
of the (unphysical) assumption, that the thickness of the 
null shell is zero, and for the collapse of the shell with 
the finite thickness this divergence becomes softer or is 
absent. To obtain a solution of the HD gravity equations 
for such a thick shell we proceed as follows [llj ■ Consider 
a set of spherical null shells collapsing to the same spatial 
point r = 0 , but passing it at different moments of time t. 
In the continuous limit, one obtains a distribution of the 
matter, that describes a spherical thick null shell which 
initially collapses and has a mass profile M(t + r), and 
after passing through the center it re-collapses with the 
mass profile M(t — r). In the linear in M approximation 
the gravitational field of such a shell can be obtained by 
averaging the solutions (EH1) . We denote by (( dh 2 )) the 
result of the averaging of the perturbation (dh 2 ). For 
simplicity we present here the expression for ((dh 2 )) for 
the case when M is constant, and the time duration of 
the thick shell is b, so that the total mass M of the shell 
is Mb. In the domain of the intersection of the incoming 
and outgoing null fluid fluxes the metric is static. The 
calculations give (see 0 for more details) 

2 GM. 9 dr 2 1, 9 . 9 . 

((dh 2 )) = --— [c 0 dt +C 2 ^- + -(cor 2 -c 2 )duj 2 ], (38) 

br r z 2 


where Ck = / . dxx k F(r 2 — x 2 ). It is easy to check that 
constant C, which enters (l27l) does not contribute to the 
curvature. For this reason we put C = 0. Using expan¬ 
sion (EH of F(z) one obtains 

c 0 = — tt[(6u - 5)S 2 + 3S], 
y 

j.5 

C2 = “ 225 [(3 ° U “ 31)52 + 15 ^’ ( 39 ) 

where u = lnr — c — In 2. For small M the function g 
remains positive, while the Kretschmann invariant is 

32 

1Z 2 ~ ^G 2 M 2 [(36u 2 + 5)Sf + 36uS 2 S + OS' 2 ]. (40) 

Hence, a collapse of a thick null shell in the theory with 
higher derivatives results in the logarithmic singularity 
of the curvature. However, if such a theory is ^-regular, 
the curvature is finite. In particular, this property is valid 
for any regular theory with higher derivatives. For the 
ghost-free theory of gravity the Kretschmann invariant 
7 Z 2 ~ ^G 2 M 2 /x 4 5 6 7 is always finite at r = 0 ( 0 . 

Let us summarize the results. Our main conclusion 
is that in the higher dimensional and ghost free theo¬ 
ries of gravity there exist a mass gap for a mini black 
hole production. Such theories possess a characteristic 
length scale /r -1 and the apparent horizon is not formed 
when the gravitational radius of the mass M, 2 GM, is 
smaller that . For a regular (and ^-regular) theory 
the Kretschmann curvature invariant for the thick null 
shell collapse is regular. Whether in the strong gravity 
limit this property remains valid, and the black hole inte¬ 
rior does not contain singularities, is an interesting open 
question. 
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